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ABSTRACT 



Context. Although a bluer-when-brighter trend is commonly observed in blazars, the opposite trend of redder-when-brighter has also 
been found in some blazars. 

Aims. We investigate the frequency of the redder-when-brighter trend in flat-spectrum radio quasars (FSRQs). 

Methods. We investigate the optical variability of 29 FSRQs in the SDSS Stripe 82 region using SDSS DR7 released multi-epoch 
data covering about nine years. We determined the spectral index by fitting a powerlaw to SDSS ugriz photometric data, and explored 
the relationship between the spectral index and source brightness. 

Results. For all FSRQs studied, we detect variations in r band flux of overall amplitude between 0.24 mag and 3.46 mag in dilferent 
sources. Fourteen of 29 FSRQs display a bluer-when-brighter trend. However, only one source exhibits a redder-when-brighter trend, 
which implies that this behavior is rare in our FSRQ sample. In this source, the thermal emission from the accretion disk may be 
responsible for the redder-when-brighter trend. 

Key words, galaxies: active - galaxies: quasars: general ~ galaxies: photometry 
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1. Introduction 

Blazars, including BL Lac objects and flat-spectrum radio 
quasars (FSRQs), are the most extreme class of active galac- 
tic nuclei (AGNs), characterized by strong and rapid variability, 
high polarization, and apparent superluminal motion. These ex- 
treme properties are generally interpreted as a consequence of 
non-thermal emission from a relativistic jet oriented close to the 
line of sight. As such, they represent a fortuitous natural labo- 
ratory with which to study the physical properties of jets, and, 
ultimately, the mechanisms of energy extraction from central su- 
permassive black holes. 

While the jet nonthermal emission is generally thought to 
be the most dominant type in FSRQs, this is not always true. 
Chen, Gu & Cao (2009a) found that thermal emission can be 
the most significant at least in optical bands for some FSRQs 
selected from the SDSS DR3 quasar catalogue. The emission of 
FSRQs is indeed a mixture of thermal emission from accretion 
disk and jet nonthermal emission. As thermal emission gener- 
ally does not vary much as clearly evident in radio-quiet AGNs 
(e.g. Ai et al. 2010), any variations in the emission of FSRQs 
could be dominated by variations in the jet nonthermal compo- 
nent. Consequently, the ratio of thermal to nonthermal emission 
varies at different source flux states. This may cause variations 
in the spectral shape where, for example, the color becomes red- 
der when the source is brighter as shown in Gu et al. (2006). In 
principle, the details of spectral variability may largely depend 



on the ratio of thermal to nonthermal emission, and the charac- 
teristics of variations in each component. 

The color/spectral behavior can be used to explore the emis- 
sion mechanisms in blazars, which remain unclear. The bluer- 
when-brighter trend (BWB) is commonly observed in blazars, 
especially BL Lac objets (e.g. Ghisellini et al. 1997; Fan et al. 
1998; Massaro et al. 1998; Ghosh et al. 2000; Clements & Carini 
2001; Raiteri et al. 2001; Villata et al. 2002; Vagnetti et al. 2003; 
Wu et al 2005, 2007; Rani et al. 2010), which can usually be ex- 
plained with shock-in-jet models. However, the opposite trend, 
i.e. redder- when-brighter trend (RWB) has also been found (e.g. 
Ghosh et al. 2000; Gu et al. 2006; Rani et al. 2010). The several 
examples are 3C 446 (Mifler 1981), PKS 0735-1-178 (Ghosh et 
al. 2000), PKS 0736-1-017 (Clements et al. 2003; Ramirez et al. 
2004), PKS 0420-014 and 3C 454.3 (Gu et al. 2006), and PKS 
1622-297 and CTA 102 (Osterman Meyer et al. 2008, 2009). 
The situation even may be more complicated because no clear 
color trend has been reported in some cases (e.g. Ghosh et al. 
2000; StaHn et al. 2006; Bottcher et al. 2007; Poon et al. 2009). 
Although there have been many studies of the color behavior of 
FSRQs, not many FSRQs have been found to display a redder- 
when-brighter trend. Moreover, while BL Lac objects display 
in general a bluer-when-brighter trend, it is unclear whether a 
redder-when-brighter trend is generally present in FSRQs (Gu 
et al. 2006; Hu et al. 2006; Rani et al. 2010). In this work, we 
investigate the spectral variability of FSRQs, to help us establish 
how common the redder- when-brighter trend is in FSRQs. 
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The layout of this paper is as follows: in Section 2, we de- 
scribe the source sample; the variability results are outlined in 
Section 3; Section 4 includes the discussion; and in the last 
section, we draw our conclusions. The cosmological parame- 
ters Hq = 70 km s"' Mpc"', = 0.3, and Qa = 0.7 ai-e 
used throughout the paper, and the spectral index is defined as 
fy oc v"'', where fy is the flux density at frequency v. 

2. Sample selection 

2.1. Quasars in Stripe 82 region 

Our initial quasar sample was selected as those quasars both in 
the SDSS DR7 quasar catalogue (Schneider et al. 2010) and 
Stripe 82 region. The SDSS DR7 quasar catalogue consists of 
105,783 spectroscopically confirmed quasars with luminosities 
brighter than M, = -22.0, with at least one emission line having 
a full width at half-maximum (FWHM) larger than 1000 km s ' 
and highly reliable redshifts. The sky coverage of the sample is 
about 9380 deg^ and the redshifts range from 0.065 to 5.46. The 
five-band (m, g, r, i, z) magnitudes have typical errors of about 
0.03 mag. The spectra cover the wavelength range from 3800 to 
9200 A with a resolution of ^ 2000 (see Schneider et al. 2010 for 
details). The Stripe 82 region, i.e. right ascension a - 20^ - 4^ 
and declination 6 = -1°.25 — i-l°.25, was repeatedly scanned 
during the SDSS-I phase (2000 - 2005) under generaUy photo- 
metric conditions, and the data are well calibrated (Lupton et al. 
2002). This region was also scanned repeatedly over the course 
of three 3-month campaigns in three successive years in 2005 - 
2007 known as the SDSS Supernova Survey (SN survey). The 
multi-epoch photometric observations therefore enable us to in- 
vestigate the optical variability of the selected quasars. 

2.2. Cross-correlation witli radio catalogues 

In this paper, we define a quasar to be a FSRQ according to 
its radio spectral index. Therefore, we cross-correlate the initial 
quasar sample with the Faint Images of the Radio Sky at Twenty 
centimeters (FIRST) 1.4-GHz radio catalogue (Becker, White 
& Helfand 1995), the Green Bank 6-cm (GB6) survey at 4.85 
GHz radio catalogue (Gregory et al. 1996), and the Parkes-MIT- 
NRAO (PMN) radio continuum survey at 4.85 GHz (Griffith & 
Wright, 1993), as well as for sources with 6 < 0°. The FIRST 
survey used the Very Large Array (VLA) to observe the sky at 
20 cm (1.4 GHz) with a beam size of 5.4 arcsec. FIRST was 
designed to cover the same region of the sky as the SDSS, and 
observed 9000 deg^ at the north Galactic cap and a smaller 2°. 5 
wide strip along the celestial equator. It is 95 per cent complete 
to 2 mJy and 80 per cent complete to the survey limit of ImJy. 
The survey contains over 800,000 unique sources, with an astro- 
metric uncertainty of < 1 arcsec. 

The GB6 survey at 4.85 GHz was executed with the 91-m 
Green Bank telescope in 1986 November and 1987 October. 
Data from both epochs were assembled into a survey covering 
the 0° <6 <15° sky down to a hmiting flux of 18 mJy, with 3.5 
arcmin resolution. GB6 contains over 75,000 sources, and has a 
positional uncertainty of about 10 arcsec at the bright end and 
about 50 arcsec for faint sources (Kimball & Ivezic 2008). The 
PMN surveys were made using the Parkes 64-m radio telescope 
at a frequency of 4850 MHz with the NRAO multibeam receiver 
mounted at the prime focus (Griffith & Wright 1993). The sur- 
veys had a spatial resolution of approximately 4' .2 FWHM and 
were made for the southern sky between declinations of -87° 
and -1-10°, and all right ascensions during June and November 



in 1990. The positional accuracy is close to 10 arcsec in each 
coordinate. The survey was divided into four declination bands. 
One of these four is the equatorial survey (-9°. 5 < 6 < +10°.0) 
covering 1.90 sr, which contains 1 1,774 sources to a flux limit of 
40 mJy and largely overlaps the GB6 survey in the declination 
range from 0° to -h10° (Griffith et al. 1995). 

The initial quasar sample was first cross-coiTelated between 
the SDSS quasar positions and the FIRST catalogue to within 2 
arcsec (see e.g. Ivezic et al. 2002; Lu et al. 2007). The result- 
ing sample of SDSS quasar positions was then cross-correlated 
with both the GB6 and PMN equatorial catalogues to within 1 ar- 
cmin (e.g. Kimball & Ivezic 2008). Owing to the different spatial 
resolutions of FIRST, GB6, and PMN, multiple FIRST counter- 
parts were found to within 1 arcmin for some quasars, although 
there is only single GB6 and/or PMN counterparts existed. The 
flux density of the closest FIRST counterparts may be dominant 
among multiple counterparts in some sources, although not in 
others. In this paper, we conservatively selected the sources with 
single FIRST counterparts to within 1 arcmin of the SDSS posi- 
tions. This selection criterion, for example, helps avoid the con- 
tamination of genuine steep-spectrum radio quasars (SSRQs). 
These SSRQs may appear as multi-FIRST components, but be 
identified as FSRQs if using the flux density of the closest FIRST 
component, which is actually much smaller than that of radio 
lobes. The optical variability of quasars with multiple FIRST 
counterparts will be presented in a forthcoming paper. 

A sample of 37 quasars was constructed from the above 
cross-coiTelations. The radio spectral index a was then calcu- 
lated between the FIRST 1.4 and either or both of the GB6 and 
PMN 4.85 GHz. Finally, 32 quasars were conventionally defined 
as FSRQs, with a < 0.5, and the remaining five quasars are 
SSRQs with a > 0.5. When counterparts were found in both the 
GB6 and PMN surveys, the spectral indices were consistent with 
each other. 

2.3. Sample 

Our final sample of 32 FSRQs is listed in Table [1] which pro- 
vides the source redshift, radio loudness, radio spectral index, 
black hole mass, disc bolometric luminosity, and the Eddington 
ratio, and the similar data of five SSRQs are listed for com- 
parison. The distribution of redshift, radio loudness, black 
hole mass, and the Eddington ratio are shown in Fig. [1] The 
source redshift is taken directly from the SDSS DR7 quasar 
catalogue, which covers 0.4 < z < 2.7. The radio loud- 
ness is from Shen et al. (2010), which ranges from log R = 
1.97 (SDSS J005205.56H-003538.1) to log R = 4.48 (SDSS 
J222646. 53+00521 1.3). However, the radio loudness was calcu- 
lated asR^ /6cm//2500, whcre /ecm and /2500 are the flux density 
at rest-frame 6 cm and 2500A, respectively (see Shen et al. 2010 
for more details). Among 32 FSRQs, 14 sources have an inverted 
spectral index between 1.4 and 4.85 GHz with a < 0. 

Black hole masses are estimated from the various empirical 
relations in the literature by using the luminosity and FWHM of 
broad H/3, Mg II, and C IV lines, i.e., Vestergaard & Peterson 
(2006) for H/3, and Kong et al. (2006) for Mg II and C IV. The 
luminosity and FWHM of broad H/3, Mg II, and C IV lines 
are adopted from the measurements in Shen et al. (2010). The 
BLR luminosity Lblr is derived following Celotti, Padovani & 
Ghisellini (1997) by scaling the strong broad emission lines H/3, 
Mg II, and C IV to the quasar template spectrum of Francis et 
al. (1991), in which Lyo' is used as a flux reference of 100. By 
adding the contribution of Ha with a value of 77, the total rela- 
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live BLR flux is 555.77, which consists of Ufi at 22, Mg II at 34 
and C IV at 63 (Francis et al. 1991; Celotti et al. 1997). From the 
BLR luminosity, we estimate the disc bolometric luminosity as 
LboI = IOLblr (Netzer 1990). The black hole mass ranges from 
to most of sources being in 10**-^ - lO'^ ^Mg 

range, which is typical of blazars (Ghisellini et al. 2010). The 
Eddington ratio log LBoi/i^Edd ranges from -1.96 to 0.57, with 
most sources being in the -1.0 - 0.0 range, which is indicative of 
a standard thin disk (Shakura & Sunyaev 1973) being present in 
all FSRQs. 

3. Results 

The SDSS DR7 CAS contains the Stripe82 database, containing 
all imaging from SDSS stripe 82 along the celestial equator at 
the southern Galactic cap. It includes a total of 303 runs, cover- 
ing any given piece of the close to 270 deg^ area approximately 
80 times. Only about one-quarter of the stripe 82 scans were ob- 
tained in photometric conditions, the remainder being taken un- 
der variable clouds and often poorer than normal seeing. For the 
runs that were non-photometric, an approximate calibration, us- 
ing the photometric frames as reference, was derived and made 
available in the CAS Stripe82 database. In this work, we directly 
use the point-spread-function magnitudes in the CAS Stripe82 
database from the photometric data obtained during the SDSS-I 
phase from data release 7 (DR7; Abazajian et al. 2009) and the 
SN survey during 2005 - 2007. The typical measurement error 
in magnitude is about 0.03 mag. 

Among 32 FSRQs, we select the sources classified as point 
sources in all observational runs. Only data with good measure- 
ments (high-quality photometry) are selected following the rec- 
ommendations in the SDSS instruction^ Moreover, we insist on 
the ugriz magnitude being brighter than the magnitude limit in 
each band, i.e. 22.0, 22.2, 22.2, 21.3, 20.5 at u, g, r, i, z, respec- 
tively. The data taken at cloudy conditions are also excluded. We 
calculate the spectral index from the linear fit to the log/v - log v 
relation after applying an extinction correction to the ugriz flux 
density and taking the flux error into account. In most cases, 
the linear fit gave good fits. Each cycle of ugriz photometry was 
usually completed within five minutes, i.e. quasi-simultaneously, 
therefore, the spectral index calculation will not be seriously in- 
fluenced by any source variations. 

3.1. Variability 

Among 32 FSRQs, three sources were excluded from our anal- 
ysis for various reasons. SDSS J015243.14-H002039.6 is classi- 
fied as a galaxy in the SDSS pipeline, although it is classified 
as a quasar in the SDSS DR7 quasar catalogue. The images of 
SDSS J215349.75+0031 19.5 are usually saturated because of a 
nearby bright object. Owing to the low quality of the u band im- 
ages, SDSS J222646.53 +0052 11.3 was also excluded because 
we aim to derive the spectral index using only high quality ugriz 
photometric data. 

All remaining 29 FSRQs show large amplitude variations 
with overall variations in r band Ar = 0.24 - 3.46 mag 
(see Table [TJ, which is much larger than that of radio quiet 
AGNs, 0.05 - 0.3 mag (e.g. Ai et al. 2010), but typical of 
blazars (e.g. Gu et al. 2006). There are four sources with 
Ar > 1 mag, i.e. SDSS JOOl 130.400+005751.80 - Ar = 
3.46 mag, SDSS J023 105.597+000843.61 - Ar = 1.02 mag, 
SDSS J0255 15.096+003740.55 - Ar = 1.70 mag, and SDSS 

' http://www.sdss.org/dr7/products/catalogs/flags.html 



J235936.817-003112.78 - Ar = 1.20 mag. In general, the varia- 
tions in different bands show similar trends, as can be seen from 
the light curve of SDSS J0255 15.096+003740.55 for example in 
Fig. |2] This source has been observed every year from 2000 to 
2007. From Fig. |2] the source remained stable in 2000 - 2001, 
gradually brightening after that. It probably reached its brightest 
state between the 2004 and 2005 observational runs, after which 
the source gradually decreased in brightness until reaching its 
faintest magnitude in 2007. The fluctuations can also be seen in 
each observational session, for example, there having been about 
0.5 mag of amplitude variations in the observational run of 2005. 

3.2. Spectral index & brightness relation 

The correlation between the spectral index ay and PSF r mag- 
nitude was checked for all sources using the Spearman rank 
correlation analysis method. We found that 15 of 29 FSRQs 
show a significant correlation at a confidence level of > 99%, 
of which 14 FSRQs show positive correlations, and only one 
FSRQ (SDSS JOOl 130.40+005751.7, see next section for de- 
tails) shows a negative correlation (see Table 1). In contrast, two 
of the five SSRQs have significant positive correlations at a con- 
fidence level of > 99% (see Table 1). 

As an example, the correlation between and r magnitude 
is shown in Fig.[3]for SDSS J0255 15.096+003740.55. A signifi- 
cant positive correlation can be clearly seen. The Spearman cor- 
relation analysis gives a correlation coefficient of r, = 0.762 at a 
confidence level of » 99.99%. This positive correlation means 
that the source spectra becomes flatter as the source bright- 
ens, in other words, the source follows a bluer-when-brighter 
trend. We noted that the spectral indices are all smaller than 
1.0, implying a rising SED in the ugriz region, more strictly, 
in 1755 - 4415 A wavelength region in the source rest-frame 
(source redshift z - 1.0228). Since the overall variation ampli- 
tude Ar = 1.70 mag is much larger than the typical values of 
radio-quiet AGNs (e.g. Ai et al. 2010), the variations are most 
likely to originate from jet nonthermal emissions. Therefore, the 
bluer-when-brighter trend implies that the synchrotron peak fre- 
quency likely increases as the source brightens. Assuming that 
the source SED in the ugriz region closely follows the overall 
SED of synchrotron emission, ay < 1.0 means that the syn- 
chrotron peak frequency may likely be Vpeak > lO''''^^ Hz (i.e. 
correspond to wavelengths shorter than 1755 A). This would in- 
dicate that this source should be classified as a high-synchrotron- 
peaked FSRQ which however is not common in FSRQs (Abdo et 
al. 2010, see also Chen et al. 2009a). Interestingly, we found that 
17 of 29 FSRQs have ay < 1 .0 all the time during variability (see 
Col. 9 of Table[T]), thus likely have a high synchrotron peak fre- 
quency if their ugriz SED closely follows the overall synchrotron 
SED. Nevertheless, the firm SED classification should be veri- 
fied by simultaneous multi-band data with a much wider SED 
coverage. 

3.3. SDSS J0011 30.40+005751. 7 

As shown in the previous section, SDSS JOOl 130.40+005751.7 
is the only FSRQ displaying a negative correlation between the 
spectral index and r magnitude. The negative correlation shows 
that the source spectrum becomes steeper when the source is 
brighter, i.e. that there is a redder-when-brighter trend. This 
source is a known FSRQ, which is included in the Radio Optical 
X-ray ASDC (ROXA) blazar sample (Turriziani, Cavazzuti & 
Giommi 2007). It appears as a compact radio source in the NVSS 
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and FIRST maps, with a flux density of /first = 156.07 mJy, and 
/nvss - 167.2 mJy in FIRST and NVSS, respectively. The spec- 
tral index calculated for FIRST 1 .4 GHz and PMN 4.85 GHz flux 
density of /pmn = 170 mJy is a = -0.07, while it is a = 0.13 
when the GB6 4.85 GHz flux density of /gb6 = 132.76 mJy is 
used. It is in the fifth VLBA calibrator survey (VCS5: Kovalev 
et al. 2007), and is compact in the 8.6 GHz VLBA map, with 
an integrated flux density of /im = 90 mJy, and a flux den- 
sity for the unresolved component of f^ms - 70 mJy. It is also 
included in the Combined Radio All-sky Targeted Eight GHz 
Survey (CRATES: Healey et al. 2007) with /8.4ghz = 278.7 mJy 
from VLA maps. 

The light curves of SDSS JOOl 130.40H-005751.7 at u, g, r, i, 
and z bands are shown in Fig. |4] The variations exhibits sim- 
ilar trends in all bands, over nine years from 1998 to 2007. 
The overall variation in r band is 3.46 mag, which is shown in 
Fig.|5] In the observational sessions, the source was brightest at 
MJD=51081 with r = 17.49, while faintest at MJD=54411 with 
r = 20.95. This source seems to vary all the time, with about 1 
mag fluctuations in most observational sessions (see Fig.jSj. 

The negative correlation between the spectral index and 
r magnitude is shown in Fig.|6] The Spearman coiTelation anal- 
ysis shows a significant negative correlation with a coiTelation 
coefficient of - -0.606 at confidence level of > 99.99%. 
This source is included in the first Fermi Large Area Telescope 
AGN catalogue with a photon index of 2.51 + 0.15 and clas- 
sified as a low-synchrotron-peaked FSRQs with Vpeak < 10'"* 
Hz (Abdo et al. 2010). With the source redshift z = 1.4934, 
SDSS ugriz wavebands coiTespond to the wavelength range of 
1424 - 3582A in the source rest-frame, which is therefore likely 
at the falling part of synchrotron SED. For a sample of 17 radio- 
quiet AGNs, Shang et al. (2005) show that the spectral break 
occurs at around 1 100 A for most objects, which is thought to be 
closely related to the big blue bump. If this spectral break also 
exists in SDSS JOOl 130.40-1-005751.7, we would expect to ob- 
serve the rising part of accretion disk thermal emission when it 
dominates over the nonthermal emission. While the spectral in- 
dex and brightness generally follows an anti-correlation, we note 
that it is more complex than a simple anti-coiTelation, as shown 
in Fig. |6] While the rising optical spectral (o-y < 10) is indeed 
present at low flux state, the declining spectra {ay > 1.0) are 
also clearly visible at low fluxes. In contrast, the optical spectra 
decline (a^ > 1 .0) at high flux states, which implies that nonther- 
mal low-peak-frequency synchrotron emission then provides the 
greatest contribution. 



4. Discussions 

The bluer- when-brighter trend is commonly observed in blazars, 
as it has been found by different investigators that the amplitude 
of the variations is systematically larger at higher frequencies, 
suggesting that the spectrum becomes steeper when the source 
brightness decreases, and flatter when it increases (Racine 1970; 
Gear, Robson & Brown 1986; GhiselHni et al. 1997; Massaro 
et al. 1998). Investigations of spectral variability have also de- 
tected this general trend (D'Amicis et al. 2002; Vagnetti et al. 
2003; Fiorucci, Ciprini & Tosti 2004). From our investigations, 
the variations in the spectral index of 14 of 29 FSRQs indeed 
follow this trend. This common phenomenon has different pos- 
sible explanations (Fiorucci et al. 2004). It may indicate the pres- 
ence of two components that contribute to the overall emission 
in the optical region, one variable (with a flatter slope Cyai, where 
fy oc V""), and the other stable (with aconst > ffvai-)- Alternatively, 



it may be interpreted using a one-component synchrotron model: 
the more intense the energy release, the higher the particle's en- 
ergy (Fiorucci et al. 2004). Another possible explanation is that 
the luminosity increase was caused by the injection of fresh elec- 
trons, with an energy distribution harder than that of the pre- 
vious, partially cooled ones (e.g. Kirk, Rieger & Mastichiadis 
1998; Mastichiadis & Kirk 2002). Finally, it could be due to 
Doppler factor variations in a spectrum slightly deviating from a 
power law, e.g. Doppler factor variations in a 'convex' spectrum 
(Villata et al. 2004). 

However, some evidence that the amplitudes of variations 
are not systematically larger at high frequencies has been found 
on several occasions (see, for example: Malkan & Moore 1986; 
Brown et al. 1989; Massaro et al. 1998; Ghosh et al. 2000; 
Clements et al. 2003; Ramirez et al. 2004). Based on their re- 
sults, Ghosh et al. (2000) suggested that it may be incorrect to 
generalize by saying that the amplitude of the variation in blazars 
is systematically larger at higher frequency. They found that 
the spectral slope of AO 0235-1-164 remained almost constant 
when its brightness increased. In S5 0716+714, while the long- 
term spectral change is evidently BWB, both BWB chromatism 
and achromatism were found in both bright and dim source flux 
states among the microvariability nights (Poon et al. 2009). The 
authors claimed that the BWB behavior can be explained by the 
shock-in-jet model, while the achromatic trend may be due to 
geometric effects. Interestingly, we did not find any significant 
correlations between spectral index and brightness for 14 FSRQs 
in our sample. 

The redder-when-brighter trend has also been found in sev- 
eral cases (e.g. Gu et al. 2006, Ramirez et al. 2004). The compe- 
tition between the thermal accretion disk and synchrotron emis- 
sion was used to explain the RWB trend in FSRQs (Gu et al. 
2006). The optical emission of FSRQs is contaminated by ther- 
mal emission from the accretion disk. Since FSRQs are usually 
low-synchrotron-peaked sources with Vpeak < 10'"* Hz (Abdo et 
al. 2010), the big blue bump of accretion disk emission will flat- 
ten the spectral slope in the optical region, resulting in a flatter 
composite spectrum than the non-thermal component. When the 
object brightens, the non-thermal component has a more domi- 
nant contribution to the total flux, and the composite spectrum 
steepens. The existence of thermal accretion disk emission is 
supported by the UV excess in 3C 454.3 during lower flux states 
(Raiteri et al. 2007). The thermal blue bump has also been ob- 
served in 3C 345 (Webb et al. 1994). Interestingly, the emer- 
gence of accretion disk emission was detected in NLSls PMN 
J0948-I-0022 as its synchrotron emission dropped (Abdo et al. 
2009). While the spectral behaviors are not homogenous across 
our FSRQs sample, it may also be inhomogenous even for a sin- 
gle source. From the four-year five-waveband monitoring, Wu 
et al. (2010) found both a bluer- when-brighter and redder- when- 
brighter trend in FSRQ 3C 345 when using a different color in- 
dex. The RWB trend can be explained by the less variable ther- 
mal components. Moreover, the authors argued that the color be- 
haviors of FSRQs are linked not only to their emission process, 
but also to the redshift and strengths of the less variable thermal 
components (Wu et al. 2010). Hence, it is most likely that the 
spectral behaviors of FSRQs depend on the position of the syn- 
chrotron peak frequency, the sampled optical wavelength range 
in the source rest frame (which depends on the redshift), and 
the positions of thermal blue bump and its strength compared to 
the jet emission. We note that FSRQs are supposed to display a 
redder- when-brighter trend, for example, in at least two of three 
sources, i.e., 3C 454.3 and PKS 0420-014 (Gu et al. 2006), and 
four of the six sources, PKS 0420-014, 4C 29.45, PKS 1510- 
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089, and 3C 454.3 (Rani et al. 2010), all of which are known 
low-synchrotron-peaked FSRQs. In contrast, our results imply 
that the redder-when-brighter trend may be rarer for FSRQs, at 
least for our present sample. However, the bluer-when-brighter 
trend is more common. 

Radio-loud quasars, particularly FSRQs, differ from radio- 
quiet AGNs in terms of their contribution of the nonthermal jet 
emission in optical continuum emission, in addition to the ther- 
mal emission from the accretion disk. In particular, jet emis- 
sion can be dominant because of the beaming effect when the 
jet is moving towards us with a small viewing angle. To in- 
spect the continuum emission in our FSRQs, we plot in Fig. |7] 
the relationship between the continuum luminosity at 3000 A 
and broad Mg 11 line luminosity for 28 FSRQs, both of which 
are adopted from Shen et al. (2010). For comparison, the ordi- 
nary least-square bisector linear fit AL^ 3000A - ^^-^ ^Ma^u ^'^^ ^ 
sample of radio-quiet AGNs in Kong et al. (2006) is also plot- 
ted with the solid line (see also Chen et al. 2009b). It can be 
seen that the 3000 A luminosity of most sources lies above the 
AL^2Qoo A ~ ^Mgii relation of radio-quiet AGNs. However, the de- 
viation of /iL^3()QQ ^ is within 0.5 dex. This excess is most likely 
to be caused by the jet nonthermal emission. The spectra of 
SDSS J001130.40-H005751.7 at MJD=51793 was used to mea- 
sure the broad emission lines and continuum (Shen et al. 2010), 
from which the corresponding r magnitude is about 19.6. If the 
luminosity at 3000 A of SDSS JOOl 130.40-H005751.7 follows 
the variation amplitude of the r magnitude, we can estimate its 
highest and lowest values, which are indicated by two rectangles 
assuming a same Mg 11 luminosity. Therefore, in the highest flux 
state, the continuum emission could be dominated completely by 
the nonthermal emission by more than one order of magnitude 
(see Fig.|7]l. However, in the lowest flux state, the thermal emis- 
sion is likely to be responsible for the continuum emission. 

While FSRQs are usually associated with core-dominated 
radio quasars, SSRQs are generally related to lobe-dominated 
ones, usually with two large-scale optically thin radio lobes. 
Usually, the beaming effect is not severe in SSRQs because of 
the relatively large viewing angle. Therefore, the optical con- 
tinuum of SSRQs could be dominated by the thermal emission. 
In Fig. I2I two SSRQs with available Mg 11 measurements all 
lie close to but below the solid line for radio quiet AGNs. This 
likely implies that the thermal emission from accretion disk is in- 
deed dominated in the optical continuum at the epoch of SDSS 
spectra. Interestingly, two of five SSRQs display a BWB trend 
(see Table[TJ- The BWB trend has also been found in radio quiet 
AGNs, and the variability was found to be anti-coiTelated with 
the Eddington ratio (Ai et al. 2010). In radio quiet AGNs, the 
optical variability could be due to the variation in the accretion 
process, for example, the variation in the accretion rate (e.g. Li 
& Cao 2008), since the jet is either weak or absent. These expla- 
nations could also be used to explain the variability of SSRQs 
if the optical emission is indeed mainly from the accretion disk. 
However, that the variation is caused by jet nonthermal emis- 
sion cannot be completely excluded in some cases because the 
variability amplitudes are generally larger that the typical values 
for radio quiet AGNs (see Table[T]). In extreme cases, the overall 
amplitude is Ar = 0.84 mag in SDSS J012401.76-H003500.9. It 
may be more likely that both thermal and nonthermal emission 
contribute to the variability. Further multi-waveband monitoring 
is needed to help resolve these uncertainties, especially the spec- 
troscopic monitoring. 



5. Summary 

We have constructed a sample of 32 FSRQs in the SDSS Stripe 
82 region. The variability and the relationship between the spec- 
tral index and brightness were investigated for 29 FSRQs. We 
found that all FSRQs show large-amplitude overall variations, 
e.g. from 0.24 to 3.46 mag at r band. We only found a sig- 
nificant negative coiTelation between the spectral index and 
r magnitude (i.e. redder-when-brighter) in one FSRQ, SDSS 
JOOl 130.40+005751.7. This implies that the RWB trend is rare 
in our FSRQ sample. This behavior could be explained by the 
contribution of thermal accretion disk emission in the optical 
wavebands. In contrast, a bluer-when-brighter trend is more 
common in FSRQs, and found in 14 of the 29 FSRQs studied 
here. 
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Fig. 1. The histogram of sample parameters: redshift (upper left), 
radio loudness (upper right), black hole mass (lower left), and 
the Eddington ratio (lower right). SDSS J222646.53+005211.3 
is not included in the histogram of black hole mass and 
Eddington ratio because of the lack of broad emission-line mea- 
surements. 



Fig. 3. The relationship between the spectral index and the PSF 
magnitude at r band for SDSS J025515.09-H003740.5. A signifi- 
cant positive correlation is present, which impUes a bluer-when- 
brighter trend. 
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Fig. 4. The ugriz band light curve of SDSS 
JOOl 130.40+005751.7 (from top to bottom). 



Fig. 2. The ugriz band light curve of SDSS 
J0255 15.09-1-003740.5 (from top to bottom). 
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Fig.5. The r band light curve of SDSS J001130.40+005751.7. 
The light curves of the second to fifth observational sessions are 
zoomed in the second to fifth panels from up to bottom. The 
dashed lines represent 1 mag variations. 




43 
42.5 



43.5 44.0 
logic L»sii (erg s"') 



Fig. 7. The plot of broad Mg II fine luminosity and continuum lu- 
minosity at 3000 A. The solid line is the OLS bisector linear fit to 
radio-quiet AGNs in Kong et al. (2006), /IL , = 78.5 LP^^^. 
The open circles are FSRQs, sohd circles are SSRQs, and open 
triangle is SDSS J001130.40-I-00575 1.7. The vertical dashed line 
represents the luminosity range at 3000 A assuming its variation 
is identical to that in the r band, which is connected by two rect- 
angles. The dotted lines show the one order of magnitude devia- 
tion of continuum luminosity at 3000 A from the solid Une. 
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Fig. 6. The relationship between the spectral index and the PSF 
magnitude at r band for SDSS J001130.40+005751.7. A signif- 
icant anti-correlation is present, which implies a redder-when- 
brighter trend. 
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Table 1. Source list: Col. 1 - SDSS source name; Col. 2 - red- 
shift; Col. 3 - radio loudness; Col. 4 - the spectral index between 
1.4 and 4.85 GHz, of which PMN flux density is used for those 
labeled with otherwise GB6 one is used; Col. 5 - black hole 
mass in unit of solar mass; Col. 6 - disc bolometric luminosity in 
unit of erg s~'; Col. 7 - the Eddington ratio / = Lbol/Ledd', Col. 
8 - overaU variation in r band; Col. 9 - the range of variation in 
the spectral index Oy-, Cols. (10-11) - the Spearman correlation 
coefficient and probabiUty level, respectively. 
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